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We calculate the quantum correction to the classical conductance of a disordered mesoscopic 
normal-superconducting (NS) junction in which the electron spatial and spin degrees of freedom 
are coupled by an appreciable spin orbit interaction. We use random matrix theory to describe 
the scattering in the normal part of the junction and consider both quasi-ballistic and diffusive 
junctions. The dependence of the junction conductance on the Schottky barrier transparency at the 
NS interface is also considered. We find that the quantum correction is sensitive to the breaking of 
spin rotation symmetry even when the junction is in a magnetic field and time reversal symmetry 
is broken. We demonstrate that this sensitivity is due to quantum interference between scattering 
processes which involve electrons and holes traversing closed loops in the same direction. We explain 
why such processes are sensitive to the spin orbit interaction but not to a magnetic field. Finally 
we consider the effect of the spin orbit interaction on the phenomenon of "reflectionless tunnelling." 
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I. INTRODUCTION 

In mesoscopic samples of disordered normal metals at low temperatures it is possible to observe a quantum correction 
to the classical Boltzmann conductance [1]. Here mesoscopic means that within the sample, the interaction of a single 
electron with other degrees of freedom, such as other electrons, phonons, magnetic impurities etc., can be neglected. 
Such a situation can be realised in semiconductor and metal nanostructures cooled to milli-Kelvin temperatures. The 
origin of the quantum correction is quantum interference between time reversed scattering paths. (By "time reversed 
path" we mean that the electron follows the same trajectory but in the opposite sense.) The correction is sensitive 
to the breaking of spin rotation symmetry and is suppressed by the breaking of time reversal symmetry [2]. The 
spin orbit interaction has the important property of breaking spin rotation symmetry while preserving time reversal 
symmetry. 

As an example, the two probe conductance G = (e 2 /h)g of a quasi-one dimensional mesoscopic wire can be expressed 
in the form g = g cl + Sg [3]. The classical conductance of the wire is O(N): g cl — N / (1 + s) where N is the number of 
scattering channels and s = L/l is the ratio of the length L of the wire to the elastic mean free path In the absence 
of an applied magnetic field the quantum correction Sg, also called the "weak localisation correction" , is O(N ). More 
exactly, in the diffusive regime 1 <C s *C N we have Sg = —2/3 if the spin orbit interaction is negligible and Sg = +1/3 
if it is not. In a magnetic field sufficient to break time reversal symmetry the quantum correction is suppressed and 
becomes 0(1/N). 

More recently, a quantum correction to the classical conductance of a disordered mesoscopic normal superconducting 
(NS) junction was observed [4]. The correction is most pronounced in junctions where the transparency T of the 
Schottky barrier at the NS interface is low V -C 1 and the length of the normal part is such that Ts ~ 1. Under these 
conditions a phenomenon known as "reflectionless tunnelling" occurs. In [5] it is suggested that quantum interference 
is responsible for this effect. The type of scattering paths involved in this quantum interference contain a segment 
where the path of an electron (hole) incident on the NS boundary is subsequently retraced by an Andreev reflected 
hole (electron.) The enhancement of the junction conductance can be as much as several times greater than the 
classical conductance. 
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In this paper we investigate the effect of the coupling of the spin and spatial degrees of freedom of the electrons by 
an appreciable spin orbit interaction has on this quantum correction. As previously mentioned, this has an important 
effect on the quantum correction in a normal metal. As the electron diffuses through the metal, the spin orbit coupling 
causes a simultaneous diffusion of the direction of its spin. Interference between different spin states is now possible 
and this modifies the interference responsible for the weak localisation correction. What role then, if any, does this 
spin diffusion play in the NS junction? 

At first sight the answer to this question would be appear to be none at all. The argument is as follows: the spin 
rotations of an electron and an Andreev reflected hole, which traverse time reversed paths, are exactly opposite and 
cancel each other. Since interference between scattering processes involving such paths is believed to be responsible 
for the quantum correction, it should be insensitive to the spin orbit interaction. As we shall demonstrate below this 
is too simple and in fact the spin orbit interaction does affect the quantum correction to the classical conductance 
of an NS junction. We have found the flaw in the argument is that it is not only processes involving electrons and 
holes traversing time reversed paths which are responsible for the quantum correction. Processes in which electrons 
and holes traverse loops in the same sense, what we shall call here identical paths, also contribute and in this case 
the spin rotations do not cancel each other. The most striking consequence of the existence of interference involving 
such paths is that the quantum correction is sensitive to the spin orbit interaction even in a magnetic field. 

In the normal metal we have seen that the quantum correction to the classical conductance is suppressed when time 
reversal symmetry is broken in an applied magnetic field. This is easily understood since the symmetry between time 
reversed paths is broken when the time reversal symmetry is broken. In the NS junction however, the contribution 
of identical paths is still present even in a magnetic field; the electron and hole carry opposite charges, and so the 
Aharonov-Bohm phases that they accumulate as they follow such trajectories cancel. Moreover, since the electron and 
hole undergo identical spin rotations, as they follow identical paths, this contribution is sensitive to the spin diffusion 
induced by the spin orbit interaction. 

The paper is organised as follows. In Sections II -IV we deal with some necessary preliminaries: the Bogolubov-de 
Gennes equations, scattering theory and conductance formulae for the NS junction, generalising the standard theory 
to take into account the spin orbit interaction. 

In Section V we investigate a quasi-ballistic NS junction, that is to say a junction whose length is shorter than the 
mean free path but long enough so that its scattering matrix is well described by a certain random phase approximation 
[6]. We express the junction conductance Gjvs as a sum of two terms 

2e 2 2e 2 

Gns = — 9ns = — (9 ns + 6 9ns) (1-1) 

a classical conductance g^ s and a quantum correction Sgws due to quantum interference. We determine the quantum 
correction for the four ensembles listed in Table I and the results are presented in Table II. In Section VI we present 
a semiclassical interpretation of these results which permits us to identify the type of scattering paths which interfere 
to produce the quantum correction. In agreement with Ref [5], we find that there is an important contribution due 
to interference between processes involving electrons and holes traversing time reversed paths, but that there is also 
a second contribution from processes involving electrons and holes traversing identical paths. We explain why this 
contribution is sensitive to the breaking of spin rotation symmetry but not to the breaking of time reversal symmetry. 

After reaching a firm understanding of the quasi-ballistic junction we consider the diffusive regime, looking first at 
junctions without a Schottky barrier at the NS interface in Section VII. In zero magnetic field the quantum correction 
is known to be of O(N ) [7-9]. In a magnetic field the authors of Ref [10] find that the correction though smaller is 
still O(N ). The effect of the spin orbit interaction was not considered in Ref [10] and so we repeat their calculation 
taking it into account. We find that breaking spin rotation symmetry multiplies the correction by a factor of minus 
one half, regardless of whether time reversal symmetry is broken or not. 

In Section VIII we discuss the dramatic enhancement of the junction conductance, known as "reflectionless tun- 
nelling," which is observable when r « 1 and Ts ~ 1. Under these conditions the contribution of time reversed 
paths is O(N) and dominates that of identical paths which is 0(N°) so that the reflectionless tunnelling effect is, in 
a first approximation, insensitive to the spin orbit interaction. This conclusion has been confirmed by carrying out a 
numerical simulation of a junction under the relevant conditions. 

II. BOGOLUBOV-DE GENNES EQUATIONS AND THE SPIN ORBIT INTERACTION 

The Bogolubov-de Gennes (BdG) equations [11] appropriate for a metal where the electrons' spatial and spin degrees 
of freedom are coupled by a significant spin-orbit interaction, are 
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where 

H e = Hq — Ep + il 

Here Hq = Hq(t, a, p) is the single electron Hamiltonian ol the metal incorporating the spin-orbit interaction and Ep 
is the Fermi energy. In deriving this equation, it is assumed that an attractive point like interaction exists between 
the electrons: 

V(f 1 ,f 2 ) = ~V 5(f 1 -f 2 ) (2.2) 

The full interacting Hamiltonian is reduced to an effective non-interacting Hamiltonian by introducing effective po- 
tentials 

A(f) = VoW^WWJ) (2-3) 



0(f,<T,/i) = Vb(l-2^)*(f, ( T)*t(F )/U ) (2.4) 

where the overline indicates a thermal average with respect to the Fermi distribution function and ^ is the usual 
field operator appearing in the second quantised formulation of the interacting electron problem. The time reversal 
operator T, which appears in the BdG equations has the form 

T = pC (2.5) 

P = i(Jy = ^ ^ Y ) ( 2 ' 6 ) 

with C the operation of complex conjugation. The eigenstates of the BdG equation describe the excitations of the 
interacting electron system. The meaning of the electron ip e and ip h wavefunctions can be seen by writing the field 
operator as 

a)=J2 o) ln - T^ h n (r, ajfrt] (2.7) 

n 

where 7„ is the annihilation operator of the excitation labelled n. The corresponding eigenvalue e„ of the BdG 
equations corresponds to the energy of the excitation. With the aid of (2.7) the occupation of the eigenstates of H e 
for a given excitation can be determined. 

III. SCATTERING THEORY FOR THE NS JUNCTION 

In this section we develop the scattering theory appropriate to the normal superconducting junction shown schemat- 
ically in Figure 1. We treat the scattering at the NS interface within Andreev's approximation. We suppose that 
any impurities in the system are in the region indicated by shading in Figure 1, and we make the approximation that 
the magnetic field is zero everywhere except in this disordered region. This is reasonable for the low fields of interest 
which affect the interference between electrons and holes in the normal part of the junction. 

A. Scattering matrices for electrons and holes 

To facilitate the explanation of the formalism, it is helpful to develop the scattering theory for a definite model. We 
have chosen a lattice tight binding model, the same model which we will use later in numerical simulations. An exact 
analogous explanation is also possible for a continuum model. 

We consider a cubic lattice and take into account nearest neighbour interactions only. We denote by ip e ( x , V, z , CT ) 
the amplitude that the electron is in an s- orbital at (x,y,z) with spin a and similarly tp h (x, y, z, cr) for a hole. We 
include in the Hamiltonian a spin orbit term which arises from the Zeeman coupling of the electron spin with the 
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effective magnetic field felt by the electron as it moves in the spatially varying potential of the lattice. We ignore the 
direct Zeeman coupling of the spin to any external magnetic field. A simple calculation shows that the Hamiltonian 
has the form: 

< x,y,z,a\H e \x,y,z,fj, > 

< x,y,z,a\H e \x - l,y,z,fi > 

< x,y,z,a\H e \x,y - 1,z,/j,> 

< x,y,z,a\H e \x,y,z- > 

where 

K ^ = VoS^ ~ V\i\o x \ a ^ 

= VqS^ - (3.2) 

An external magnetic field B, applied in the +y direction is modelled by Peierl's factors in the matrix elements 
between nearest neighbours in the z direction. If a is lattice constant a = 2-kBci 2 /</> where <p = h/e is the flux 
quantum. The transverse dimensions are 1 < x < L x and 1 < y < L y . The Hamiltonian (3.1) can be regarded as a 
three dimensional generalisation of that proposed in Ref. [12] as a model for a two dimensional electron gas formed 
at the surface of a III-V semiconductor. 

The relevant energy scale of the model is determined by Vq + V 2 . For convenience we shall set this to unity with 
the choice 

V = cos0 Vi = sin (9 (3.3) 

Varying the angle 9, we may set an arbitrary ratio of normal potential coupling Vo to spin orbit coupling V\, while 
keeping the extent in energy of the density of states roughly constant. With this choice, v x ,v v and v z are all elements 
of SU(2). Using the homomorphism of SU{2) with the three dimensional rotation group SO(3), we can interpret 
the i>'s as rotations of the spin of the electron as it moves between nearest neighbours [13]. A product of nearest 
neighbour matrix elements along a path will have the form exp(i$)v where $ is the Aharonov Bohm phase picked up 
by the electron, and v € SU (2) is the rotation of the electron's spin as it traverses the path. 

The Hamiltonian H e has the form given in (3.1) everywhere except in the disordered region located in < z < L. 
There, some or all of the Hamiltonian matrix elements are supposed random. In Ref. [12] the diagonal elements were 
assumed to be independently and identically distributed, while the parameter 9 controlling the spin orbit interaction 
was held fixed. In Refs. [14,15] a random spin orbit interaction was also considered. For the present purpose we do 
not need to specify the precise distribution. 

First, we consider the scattering of electrons incident at an energy E = E F + e. To the left of the disordered section 
we expand the electron wavefunction ip e in terms of the Bloch states of (3.1) with energy E. 

ip e (x,y,z,a) = ^2 a +n ip +n (x,y,<r)cxp(+ik n z) + ^ a- n if>- n (x, y, a) ex.p(-ik n z) (3.4) 

n;lmk n <0 n\lmk n — 



V. 



cxp(-iax) v z a41 



y'" (3-1) 



tp- n (x, y, cr) = P<y^fp+n(L x -x + l,y,fi) (3.5) 

As z — > — oo, far from the disordered region, we impose the boundary condition that the allowed states consist 
exclusively of incoming and outgoing propagating waves. Thus in z < states with Imk n > are excluded. We 
denote by 2N the number of "open channels", i.e. states with lmk n = that carry a positive probability current 
in the +z direction; there are an equal number carrying current in the — z direction. We label these states so that 
+n carries a current in the +z direction and — n a current in the — z direction. After a suitable normalisation of the 
transverse wavefunctions (see Appendix A) the electric current due to electrons at the left of the disordered section is 

I e = Z j^ W+n\ 2 - \a-n\ 2 (3.6) 

n\lmk n — 

The boundary condition as z — > — oo, imposed above, ensures that only open channels, and not "closed channels" with 
Imfc„ ^ 0, contribute to the current. A similar expansion may be made on the right in terms of a set of coefficients 
{a' +n ,a'_ n } with the boundary condition that we admit only those states with lmk n > 0. Thus, far to the right of the 
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disordered section, as z — > +00, the allowed states again consist exclusively of incoming and outgoing propagating 
waves. 

The AN x AN scattering matrix for electrons S e relates the AN incoming flux amplitudes at the left, a + = 
{a + „;Imfc„ = 0} and the right a'_ = {a'_ n ;Imk n = 0} with the AN outgoing flux amplitudes at the left 
a_ = {a_„; Imk n = 0} and the right a' + = {a' +n ; Imk n = 0} 



a'_ 



The matrix S e has the structure 



t e r' 



(3.7) 



(3.8) 



in terms of the 2N x 2Areflection and transmission matrices for left incidence (r e ,t e ) and right incidence (r' e ,t' e ). 

Since we are considering time independent scattering, the currents to the left and right of the disordered section 
must be equal, and therefore it follows that S e is unitary. There is an additional restriction on S e when, in the absence 
of an applied magnetic field, the Hamiltonian is time reversal invariant ie. [H e ,T] = with T given in Eq. (2.5). For 
a suitable choice of transverse wavefunctions (see Appendix A) Se will then satisfy 



S = 



I2JV 






-12JV 



-1 2 jv 
1 2 n 



where I2N means the 2N x 2N unit matrix. This can be written in the equivalent form 



r 
' e 

(t' e ) T 



(3.9) 



(3.10) 



The simplicity of these relations, compared with for example those of Ref [16], is related to the presence of p in Eq. 
(3.5) (see Appendices A and B) 

We now turn to the scattering matrix St for the holes. From the BdG equations we can see that the hole wave- 
function ip h evolves according to 



where Hh is given by 



H h [+B] = TH e [+B]T = —H e [—B] 



(3.11) 



(3.12) 



If ij} h describes the scattering of a hole with excitation energy +e in a field +B then Tip h describes the scattering of 
an electron at energy — e also in field +B. We shall make use of this in two ways. 

Firstly, outside the disordered region B = and Hh = —H e . Thus outside the disordered region it is useful to 
expand ip h in terms of the Bloch states of H e at energy Ep — £■ At the left, for example 



ip h (x,y,z,a) = 

n;Imfe fl <0 



b- n ip +n {x, y, a) exp(+ik n z) 



E 

n;lmk n —0 



b +n ip- n (x, y, a) exp(-ifc„z) 



(3.13) 



Note that since in this region Hh = —H e the probability currents are reversed by comparison with the electron 
case. We therefore associate the coefficient the flux amplitude for a positive hole probability current in the +z 
direction, with the wavefunction proportional to exp(— ik n z). The holes carry an opposite electric charge to that of 
the electrons so that they carry an electric current 



h=- R J2 \b+n\ 2 -\b-n\ 2 (3.14) 
n;lmk n =0 

By definition, the AN x AN matrix Sh relates incoming hole probability currents to outgoing probability currents 

6_ 



S h 



b+ 
b' 



(3.15) 
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The matrix Sh has a structure similar to that of Eq. (3.8); i.e. 



S h 



rh t' h 
th r' h 



(3.16) 



Secondly by rewriting Ti[> h in terms of electron flux amplitudes, and recalling that these amplitudes are related by 
5 e , we arrive at a relation between S e (— e, +B) and Sh(+e, +B) 



S h (+e,+B) 



1-2N 







-12JV 



or 



r h (+e,+B) 
r' h (+e,+B) 
t h {+e,+B) 
t' h (+e,+B) 



S* e (-e,+B) 



[r e (-e,+B)}* 
K(-e,+B)}* 
[t e (-e,+B)}* 



-Un 
Un 



(3.17) 



(3.18) 



Again we assume here a suitable choice of transverse wavefunctions (see Appendix A.) Following a similar line of 
argument it is possible to demonstrate a relationship between S e (+e, +B) and S e (+e, —B). 



B. Andreev scattering at the NS interface 

In this section we outline the calculation of the coefficients of Andreev reflection [17] at the NS interface. These are 
essentially unchanged by the introduction of a spin orbit interaction in the materials forming the junction. We assume 
throughout that Ao <C Ep, a condition which is realised in practice. 

Far from the NS interface in the normal metal the superconducting gap A — > 0. On the other hand, far from 
the NS interface in the superconductor, A — ► A cxp(i(^) where A is real. In general the reflection coefficients will 
depend on the precise form of A in the transition region near the junction. For a point contact junction, however, it 
is permissible to assume a simple step model 



A = z(j, S 

A = A exp(iip) z e S 



(3.19) 



We are interested in the situation where the energy E = Ep + e of the incident electron is in the energy gap of 
the superconductor Ep < E < Ep + A . Anticipating somewhat in order to avoid unnecessary algebra, we find the 
electron is mainly reflected as a hole like excitation. A solution of the BdG equation in the normal metal corresponding 
to this is 







' 1 " 










exp(ik^-h)^(x, y, a) + 







he 



exp(ifc£ 



- ] z)ip { +1 l(x,y,(j) 



(3.20) 



where r£ e denotes the matrix of Andreev reflection amplitudes and the superscript (±) refer to Bloch states with 
energies E^ — Ep ± e. The first term describes an excitation where an electron above the Fermi level is incident 
from the left in channel n. The second term corresponds to an excitation in which an electron below the Fermi level is 
annihilated, i.e. to a reflected "hole" with opposite velocity and spin to that of the incoming electron. Since A <C Ep 
we can to a good approximation ignore the difference between and k„ ' and similarly the differences between 
the transverse wavefunctions. Requiring that the wavefunction and its derivative be continuous at the boundary of 
the superconductor leads to 



' he 



= i exp(— iip) , 



' eh 



(3.21) 



Here we have made the further assumption that e <C Ao, the limit of interest in what follows, and we have also given 
the reflection coefficient for an incident hole. 



G 



C. The Scottky barrier at the NS interface 



In real NS junctions, a mismatch between the conduction bands of the two materials which make up the junction 
results in the creation of a Schottky barrier at the interface. This barrier plays an important role in the physics of the 
device and so we must take it into account. We shall model the Schottky barrier as a simple planar potential barrier. 
At the barrier an incident particle may be either transmitted without a change of momentum or specularly reflected. 
We neglect any dependence of the reflection and transmission probabilities on the momentum of the incident particle 
so that the properties of the barrier are described by a single parameter re [0,1], its transparency. The transmission 
and reflection matrices which make up the electron scattering matrix Sf of the barrier are 

tf = Vr i 2N 

rf = VT^T w B (6 - ZZ) 



r' r B = -vT-l W f 



The precise form of the 27V x 27V matrix wb depends on the choice of transverse wavefunctions. For subsequent 
analysis we need only note, however, that wb is antisymmetric and unitary. The hole scattering matrix Sff for the 
barrier is related to Sf in the usual way by (3.18). 



D. Combination of scattering matrices 



It is the purpose of this section, having considered above the scattering matrices for the various components of the NS 
junction, to explain how the scattering matrices may be combined to find the total scattering matrix. We consider 
first a junction without a Schottky barrier. For the normal part we can write 



r 


t' ' 




c+ 




c_ 


t 


r' 




c'_ 




A. 



where 



are 47V-dimensional matrices and 



r e 
r h 



etc. 



c± = 



a± 




"<4" 






. b '±. 



(3.23) 



(3.24) 



(3.25) 



4iV-dimensional vectors, in the notation of Sec. IV. For the Andreev part we define the 4N x 4N scattering matrix 
S A by 



S A c' + 



By reference to Section III. this has the form 

S A = 





he 1 ™ 



rtl 







For the combined system, the 47V x 47V S matrix is defined by 

Sc + = c_ 



(3.26) 



(3.27) 



(3.28) 



In equation (3.23) we replace c'_ in terms of c' + as in (3.26); we then perform the matrix multiplication and, from the 
two resulting equations, eliminate c' + . From the expression relating c + and c_ we extract the 4iV x 47V 5* matrix of 
(3.28) to obtain 



t'S 



1 

1 - r'S A 



t 



(3.29) 
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In the electron-hole spaces, S of Eq. (3.29) has the structure 



S = 



r ee Teh 
The T hh 



(3.30) 



r ee etc. being 2N x 2N matrices. To determine the conductance we shall need the submatrix r^ e . From Eqs. (3.29) 
and (3.30) we find 



The = t' h [S 



1 



1 - r'S A 



he 



Using the structure (3.27) of S A we can write r he as 



r he — t'h r he 



1 



l-r'S A _ 

For any nonsingular operator D we now use the operator identity [18] 

1 \ 



D ee ~ D e hJJ^D he / 



to find 



rhe — t' h r he 



[l r 'e r he r 'h T th\ te 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



For a junction with a Schottky barrier we first consider the composition of the Schottky barrier and the supercon- 
ductor. The scattering matrix for this system can be obtained from Eq. (3.29) where r, t and t' are taken from the 
model (3.22) for the barrier. With the aid of the identity (3.33) and the following one [18] 



we find 



he 



D hh 



D he 



1 



Dee ~ Dehjh-D he 



r BS r BS 
' ee 'eh 

'he ' hh 



(3.35) 



(3.36) 



r BS _ 
r he 



= iexp(-i<p) [r/(2-r)] i 2N 



r BS _ 
r eh — 
r BS 

ee 
r BS 
r hh 



-(r[ B e S] r 

2 [>/r=r/(2 - r)] w B 

-(rl B e S] y 



(3.37) 



The scattering matrix for the complete system of normal part, Schottky barrier and superconductor can now be 
obtained from Eq. (3.29), where r', etc. refer to the normal metal as before, but S A is replaced by S BS of (3.37). 
(Note that in deriving Eq. (3.29) the structure (3.27) of S A was not used.) 



IV. CONDUCTANCE FORMULAE 



We assume that the bias voltage is small in comparison to A /e, and that the size of the superconducting part is 
long enough so that there is no quasi-particle current in the superconductor. In this case, the zero temperature dc 
conductance Gns of the normal-superconducting junction can be described by a simplified "Landauer" formula which 
has been derived in Ref. [19-21] 



2e 2 2e 2 
Gns = —9ns = -^-trr he rj[ e 



(4.1) 



where r^ e is the 2N x 2N matrix of electron-hole reflection amplitudes for the composite system. 
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There is an important simplification if the Hamiltonian of the system is time reversal invariant and the bias voltage 
is small in comparison to the Thouless energy E c [22], so that the energy dependence of S e can be neglected. The 
conductance then has the form [23] 



Gns 



h 



2JV r 

£ 

n=l 



T 



l 2 



2-T„ 



(4.2) 



where T n are the eigenvalues of t e t\. We have verified that this result still holds when spin rotation symmetry is 
broken by the spin orbit interaction. 

In what follows we wish to compare the quantum conductance of the NS junction, calculated from (4.1), with 
the classical conductance of the junction. This latter quantity is determined using the classical rule of combining 
conductances 



1/3 - + 1/52 



(4.3) 



This corresponds to the addition of flux intensities as opposed to flux amplitudes. The conductance associated with 
the electron traversing the normal part is 2N/s and similarly for the hole in the traversing the normal part in the 
opposite direction. The conductance of the barrier is 2N\r^f\ 2 , so that the classical conductance is 



9ns 



2iV|rf s ' 2 



he 



1 



2s\rt S V 



(4.4) 



The classical conductance is insensitive to the breaking of time reversal and spin rotation symmetries. 



V. QUASI-BALLISTIC JUNCTION 

In this section we shall calculate the conductance of an NS junction to first order in s = L/l, where L is the length of 
the normal part of the junction and I is the elastic mean free path. We neglect terms of order s 2 and above so that 
result is strictly applicable only in the limit that s« 1. Nevertheless we shall see that the results of the calculation 
shed considerable light on the origin of the quantum interference in the device. 

The 2N x 2N reflection matrix r^ e for the system consisting of the normal metal, barrier and superconductor can be 
obtained as discussed in Section HID. After expanding to second order in r' e , r' h , which is sufficient for an evaluation 
of Qns to first order in s, we have 

„. . f i r BS f , f l r BS r l r BS+ ,4.1 r BS r l r BS+ , 4.1 r BS „/ r BS _/ r BS f 

I he — he e h he e ee e ~r Ifr' frfi ' h he e ' h he e ee e ee e /r i\ 

, 4.1 r BS r i r BS r i r BS f , f i r BS r i r BS r i r BS f , f i r BS r i r BS r i r BS f , 

' L h' hh ' h' he 'e'ee L e T ^h' hh ' h' hh ' h' he Le ^ l h' he ' e< eh ' h' he te ' ■ ■ ■ 

The conductance is found by substituting this into (4.1) and performing an average over an ensemble of scattering 
matrices S e , describing a set of microscopically different but macroscopically equivalent configurations of impurities 
in the normal part of the junction. 

In principle, the distribution for S e should be calculated from some model distribution of Hamiltonians. We shall 
not, however, attempt to do that here. Instead we will assume that the resulting ensemble of scattering matrices S e 
is distributed according to the "local maximum entropy model" [3,24,25]. If the geometry is quasi-ld and the number 
of channels N sufficiently large, then results obtained with the aid of the local maximum entropy model are known 
to be identical to those obtained from the class of microscopic models described by the the nonlinear sigma model 
[26,27]. 

The distribution of S e in the local model depends on TV, s and the symmetry of the Hamiltonian, i.e. whether or 
not time reversal symmetry is broken and whether or not spin rotation symmetry is broken. There are four ensembles 
(Table. I.) The critical strengths of the magnetic field and the spin orbit interaction separating the various ensembles 
should be similar to those associated with the weak localisation effect in normal metals. The details of the distribution 
of S e for the four ensembles can be found in Appendix C. 

It will be helpful, when we come to discuss the interpretation of the results, to write (5.1) in the form 

oo 

rhe = ^2pi (5.2) 

i=l 

Each term in the series represents the contribution of a particular scattering processes to n ie . The classical conductance 
is obtained by ignoring interference between different processes and summing intensities 
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9ns =te\Y, PiP i 



(5.3) 



The quantum correction to this classical conductance is found by summing the interference between different processes 



Sg 



NS 




After carrying out the average we find that 



9ns 



2N\r^(l-2\r^\ 2 s + 0(s 2 )) 



(5.4) 



(5.5) 



which agrees with the expansion of (4.4) to the order we are considering. The explicit expressions for Sqns are 
collected together in Table II. The function / of the barrier transparency, which appears in the table, has the explicit 
form: 



/(T) 



2r 2 



(2-r) 2 



2r 2 



(2-r) 2 



i 



(5.6) 



There are two obvious limiting cases: r = 1 corresponding to a junction without a Schottky barrier and r « 1 
corresponding to a junction with a high Schottky barrier. The first point to note is that the quantum correction is 
of O(N), the same order as the classical conductance, in zero field. Secondly for T <C 1, the conductance increases as 
disorder is added to the junction. This is the essence of the dramatic reflcctionlcss tunnelling effect which we discuss 
in Section VIII. Thirdly when time reversal symmetry is broken by the application of a magnetic field the quantum 
correction is 0(N°) and not 0(1/N) as might have been expected by analogy with the weak localisation effect in 
a normal metal. The final, and perhaps the most surprising result, is that in a magnetic field the breaking of spin 
rotation symmetry by the spin orbit interaction multiplies the quantum correction by a factor of minus one half even 
though the symmetry of the Hamiltonian, in the sense of random matrix theory, is unchanged and remains unitary. 



VI. SEMICLASSICAL INTERPRETATION 



The importance of quantum interference between processes in which the path of an electron (hole) incident on the NS 
boundary is subsequently retraced by an Andreev reflected hole (electron) was first pointed out in [5] . In the absence 
of a magnetic field, and if the bias voltage is small enough, electrons which move along a path in one given sense 
are phase conjugated with holes traversing the time reversed path. In a magnetic field, or if the bias voltage is large 
enough, this phase conjugation is destroyed. However we have seen that there is a significant quantum correction even 
in a magnetic field. There must therefore be an additional source of quantum interference which is not sensitive to the 
breaking of time reversal symmetry. As we shall see the relevant processes involve paths in which an electron (hole) 
and an Andreev reflected hole (electron) traverse a loop in the same sense. In order to remain concise, we shall refer to 
such processes as containing identical paths. The interference involving such paths can only be destroyed by applying 
a large enough bias voltage. The physical importance of the bias voltage as a "symmetry breaking parameter" is 
discussed in [10]. 

Only some of the terms in (5.4) are found to be nonzero after averaging, so that in fact 

SgNS = tr (pipl + pip\ + p 5 p\ +p 7 p\^ + 0(s 2 ) (6.1) 

where pj means the jth term in (5.1.) Consider first the interference between process p\ and p§. An example of a 
scattering path which contributes to process p\ 

Pi = t' h r^t e (6.2) 

is illustrated in Figure 2. Since we are working to first order in s it is sufficient to consider the motion of the electron 
and holes, traversing from one side of the normal part of the junction to the other, as ballistic, so these trajectories 
appear as straight lines in Figure 2. Examples of paths which contribute to p$ 

P5 - t' h r^r' h r^r' e rl s t e (6.3) 

are illustrated in Figures 3 and 4. To first order in s the interference between processes p\ and P5 is 
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tr ( PlP i + p 5 p\) = 2|r£f | 2 |rf/| 2 tr (r' h r' e ) (6.4) 

The remaining terms involving pi and p-j contribute 

tv(p lP \+p 7 p\) = -2|rf e s | 4 tr(r^> (6.5) 

In the interest of brevity we will concentrate on the interference between p\ and p§. A very similar analysis is 
possible for the interference between processes p\ and pr leading to identical conclusions. We proceed by relating the 
product of electron and hole reflection matrices in (6.4) to a product of an electron and a hole Greens function. To 
simplify the algebra we shall suppose that both the spin orbit interaction and the magnetic field are zero everywhere 
except in the disordered region. We shall also suppose that L y = 1 and impose periodic boundary conditions in the 
x direction. The Bloch states at energy E have the form 



ip2m(x, z,a) = e-Kp(ik% m x)ex.p(ik 2m z)5 aiCr2m a 2m = T 

^2m+i{x,z,a) = exp(i/cf m+1 x)exp(ifc 2m+ i2;)(5 (T;(J2m+1 <? 2m +i = I 

i>-2m(x,Z,<T) = -CXp(-ik% m x)exp(-ik 2m z)5v,a- 2m <T-2m = I 

i>-(2m+i){x,z,a) = exp(-ifcf m+1 x)cxp(-ifc 2m+1 z)5 (T!CT _ (2m+1) a_ {2m +i) = T 



(6.6) 



where 



k 2m = fcfm+i = y~ m = 0, . . . , L x — 1 (6.7) 

Ltr, 



and the energy and the momenta are related by 

E = 2 cos k x m + 2cosk m (6.8) 

The reflection matrices for electrons and holes can be related to the corresponding Green's functions as indicated in 
Appendix A. 



-i\J A: sin k m sinfc„ exp(+i(fc m + k n )L) 

^*+ m ( x b^)Gt( x b,L,c7;x a ,L,a')tp^ n (x a c7') (6.9) 



x a ax b a' 



[r' h )nm = -i\/4: sin k m sin k n cxp(-i(k m + k n )L) 

X! i!>*- n {xbo)Gj l {x b ,L,(j]x a ,L,(j')il> +m (x a (j') (6.10) 

Note that, for convenience, the states (6.6) have not been normalised to carry identical currents, proper account of 
this has been taken in the expressions (6.9) and (6.10) for the reflection matrices. 

The quantum correction involves a trace over the product of r e and rv Using the relation with the Green's function 
this can be separated into two elements: an integration over the cross section involving the transverse wavefunctions 
and an average of the product of an electron and a hole Greens function. We will consider the second element first. 
This involves evaluating 

^a^ n ,a +m (x a ,X b ) = (G% (x a , L, (7_„ ; X b , L, Cr +m )G+ (x b , L, U +m ] X a , L, fT_„)) (6.11) 

Within the semiclassical approximation, as is explained in Ref. [28], Ch. 12 and Ch. 13, we can express the Greens 
functions as summations over paths: 

%b,L,<T +m ;x a ,L,(T- n ) — ^2j :Xa ^ Xb Aj exp(iSj + i^j)[vj]cr +m ,"-n {rto\ 
Gl(x b ,L,cj- n ;x a ,L,a +m ) = - T,j-. Xa ^ Xb A j exp(-iSj - i$j)[vj]a_ n ,a +m 

Substituting (6.12) into (6.11) and taking the disorder average we find that only two contributions remain: 
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A. Time reversed paths 



The path of the electron moving between x a and Xb is retraced by the Andreev reflected hole as illustrated in Figure 
3. The electron and hole charges are opposite but they traverse the path in opposite directions so the Aharonov-Bohm 
phase factors they accumulate do not cancel each other out. Thus this term is important only when the magnetic 
field is negligible and the Hamiltonian is time reversal symmetric. 

Let us consider first the situation where the spin orbit interaction is negligible. In the quasi-ld limit, the average 
of the product of the hole and electron Green's functions is independent of x a and Xb so that 

K- n ,* +m (x a ,x b ) = -\5 a _ n , a+m (6.13) 

After integrating the product of transverse wavefunctions over the cross section we find 

{[r' h ]mn[r' e ]nm) = 4 sin k m sin k n \5 a _ n ^ +m + LP. (6.14) 

where LP. denotes the contribution due to identical paths which we will consider below. We determine A by demanding 
consistency with the local maximum entropy model. In this model S e is distributed isotropically (see Appendix C) 
so that averages, such as that in (6.14), are independent of the channel indices n and m. The extent to which the 
isotropy assumption is reasonable for microscopic models is discussed in Rcfs. [29,30]. In the present case, to maintain 
consistency with the results of Section V, we are required to make the approximation that the group velocities 
appearing in (6.14) can be replaced by an average velocity. This is a reasonable assumption if the energy is such that 
we are far from any subband edges, but is of questionable validity for energies close to a subband edge. Having made 
this approximation and determined A by comparison with Section V we find 

< {r' h )2n+l,2m(r' e )2 m ,2n+l > = < (r' h hn,2 m +l(r' e )2m+l,2n >= N + 1 (6.15) 

If the spin orbit interaction is appreciable we assume that the spin direction of an electron or hole traversing a path 
is completely randomised. The matrix v describing this rotation can then be assumed to be uniformly distributed, 
relative to the invariant measure, on the group SU(2). The result is 

< (r' h ) n , m (r' e ) m ,n >= (6.16) 

After summation over all channels the contribution of this term to the quantum correction is essentially independent 
of the strength of the spin orbit interaction. This is because the spin rotation experienced by an electron traversing 
a path in a given sense cancels that of the hole traversing the same path in the opposite sense. To demonstrate 
this explicitly it suffices to note that after summation over spin indices the contribution of time reversed paths is 
proportional to trv^v — trl 2 , which is a constant independent of the distribution of v. 



B. Identical paths x a = x b 

The electron path in the normal metal contains a loop which is retraced in the same sense by the Andreev reflected 
hole as illustrated in Figure 4. The Aharonov-Bohm phase accumulated by the electron and hole, are now exactly 
opposite and cancel. Therefore this term is important even if time reversal symmetry is broken. 
Proceeding in a similar way to that above we find that the contribution of identical paths is 

(KlnmKlmn) = 4 SHI fc TO SHI fc„A5 TO ,„ (Va_ n ,cr +m ^a +m , CT _ n ) + T.R.P. (6-17) 

(Here T.R.P. refers to the contribution from time reversed paths.) The presence of the delta function, arising from the 
integration of the transverse functions over the cross section, means that identical paths only contribute for electrons 
and holes which are back scattered by the disordered region as is indicated in Figure 4. Note the novel feature that 
this is true even if time reversal symmetry is broken in the disordered region. The average over the matrices v is 
calculated with respect to the distribution p(v) = 5(v — I2) if the spin orbit interaction is negligible, and p(v) uniform 
on SU(2), if it is appreciable. As above we determine A by reference to Section V. If the spin orbit interaction is 
negligible we find: 

< ( r 'h)2n+1.2m(re)2m.2n+l > = < K) 2n 

,2m+l(Pe)2m+l,2n >~ ^m.n M : (6.18) 
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in zero magnetic field and: 

< (r' h )2n+l,2rn(r / e )2m,2n+l > = < {r' h )2n,2 m +l(r' e hm+l,2n >= ^m.nj^ (6.19) 

in a nonzero field. If the spin orbit interaction is appreciable we find 

< {r' h )n,m{r' e ) m ,n >= W^TT (6.20) 



in zero magnetic field and 



< {r' h )n,m( r 'e)m,n >= "W^y (6.21) 



in a nonzero field. The change of sign when spin rotation symmetry is broken comes about because an electron and a 
hole traversing the same path in the same direction undergo the same spin rotation. After summation over the spin 
indices, the contribution from identical paths is thus proportional to trw 2 which is sensitive to the distribution of v 

trv 2 = +2 p(v) = S(v - 1 2 ) 

= -1 p(v) uniform on SU(2) [ ° > 



VII. DIFFUSIVE JUNCTION WITHOUT A SCHOTTKY BARRIER (r = 1) 

For convenience we will defer discussion of a junction for which the Scottky barrier has a low transparency (r <C 1) 
and which exhibit reflectionless tunnelling, to Section VIII, and focus in this section on junctions without a Schottky 
barrier. 

We now have to sum the perturbation series (5.1) to all orders in s to find the quantum correction for a diffusive 
junction. In the absence of a magnetic field, this amounts to the evaluation of the integral: 



9ns = I ATp{T) 
Jo 



2-T 



(7.1) 



Here p(T) is the density of the eigenvalues of the matrix t e t\, where t e refers to the composite system composed of 
the normal metal and Schottky barrier. 

Within the local maximum entropy model, the development of the density p(T) as a function of s is described by 
a nonlinear diffusion equation [31]. The initial condition on p(T) at s = is: 

p s=0 (T) = NS(T - 1) (7.2) 

In the diffusive regime the solution to this equation is independent of the ensemble up to corrections of O(N ) [32]. 
Evaluating (7.1) we find an identical O(N) contribution to g^s for the Orthogonal and Symplectic Ensembles, while 
the following O(N ) term is ensemble dependent. For the details, we refer the reader to [32]. With the aid of this 
reference we obtain the results for the Orthogonal and Symplectic Ensembles given in Table III. 

For a diffusive junction in a magnetic field the summation of the perturbation series requires the performance of a 
number of averages over the unitary group. The authors of Ref. [10] find for the Unitary I Ensemble and T = 1 that 

27V 8s 3 + 12s 2 + 12s 
9NS = TT2s" 3(1 + 2s) 3 + ^ 

This gives the same result as Section V, in the limit that s -c 1 and yields Sgws — — 1/3 in the diffusive regime 
(1 <SC s <C N.) As for the quasi-ballistic junction, the quantum correction is 0(N°), not 0(1 /N), in a magnetic field. 

The analysis of Ref [10] can be extended to the Unitary II Ensemble. In the absence of a Schottky barrier we have 
from (3.34)' and (3.18): 

OO 

r he - i expH^i'* l-<r' e *] P *e (7.4) 
Then with the aid of (CI) and (4.1) we obtain for the conductance: 
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9ns = tr 



OO OO 

[-uVRu*VR P uTu1J2 [-\/fiu T V# 



. p=0 q=0 



(7.5) 



where u = u\u2 is uniformly distributed on U(2N). Following Ref [10] the ensemble average is accomplished in two 
steps. First the average on the unitary group and then the average over the distribution of T. The final result for the 
Unitary II Ensemble and T = 1 is: 

27V 4s 3 + 6s 2 + 6s ^ ,„„ N 
^ = TT^ + 3(1 + 2 S )3 ^ 

which reproduces the result of Section V for s«l and yields Sg^s = +1/6 in the diffusive regime. 

The most interesting result is that, like the quasi-ballistic junction, the quantum correction is sensitive to the spin 
orbit interaction even if time reversal symmetry is broken. Referring to Section VI, it seems reasonable to assume 
that it is interference involving identical paths which is responsible for this sensitivity. 

In order to confirm that quantum correction is sensitive to the spin orbit interaction, even if the junction is a 
magnetic field we have a carried out a numerical simulation. As a model for the junction we use the Hamiltonian of 
Section III. The Fermi energy Ep is measured from centre E of the energy band of the model. Since the zero of 
energy is arbitrary we can set Eq = for convenience. In the normal part of the junction, < z < L — 1, we take the 
diagonal Hamiltonian elements as: 

< x,y,z,a\H\x,y,z,n > = U(x,y, z)S <r , li < z < L z - 1 

(7.7) 

<x,y,z,a\H\x,y,z,fj,> = U 5 a ^ z = L 

The potential U(x,y,z) is the random potential due to impurities and Uo the height of the Schottky barrier. For 
simplicity we shall assume a spatially uncorrelated potential with a distribution P(U) = p(U)dU 

P (U) = 1/W -W/2 <U< +W/2 

(7.8) 

p(U) = 0, otherwise. 

To the right, in z > L + 1, is the superconductor. The conductance is calculated using the standard Green's function 
iteration technique to produce S e and then using the formulae of Sections III and IV. 

In Figure 5 we present the data for the junction conductance versus the parameter 0, defined in Section III, 
controlling the spin orbit interaction. The transverse dimensions of the junction are L x = L y = 5 and the longitudinal 
dimension L = 25. The Fermi energy Ep = 0. For these parameters N = 17. We are interested in the limit T = 1 
so we set Uq = 0. A magnetic field corresponding to a flux of l/25^o per lattice cell is applied perpendicular to the 
junction corresponding to five flux quanta through the device. The distribution of the random potential has width 
W = 3 corresponding to an estimated mean free path I ~ 4.3. The junction is in the diffusive (L z ^> I) quasi-l(i 
(L x ,L y < I, L x ,L y -c L z ) regime. The spin rotations between lattice sites are not random but differ depending on 
the direction of the displacement. Since these rotations do not commute the spin direction of the electron is effectively 
randomised and so, for 9 large enough, the Symplcctic Ensemble should be appropriate. 

As the magnitude of the spin orbit interaction is increased the junction conductance Gns increases, while the 
normal conductance Gn remains constant. This is exactly the behaviour predicted. The value SGns — 0.8e 2 /h 
obtained in the simulation is also in reasonable agreement with theoretical value 5Gns — e 2 /h. 

VIII. REFLECTIONLESS TUNNELLING 

This is an effect observed in diffusive junctions when r< 1. In the absence of a magnetic field we may calculate the 
conductance in the same way as for the diffusive junction with T = 1 but with the initial condition: 

p s=0 (T) = NS(T-r) (8.1) 

instead of (7.2). At O(N) the result is the same for both the Orthogonal and Symplectic Ensembles: 

9ns — NT 2 /2 + NsT 2 1 « s « NTs « 1 

9NS = n/( s + i/r) i < s < A^ ; r s > i ^) 
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We expect ensemble dependent corrections to this result at O(N ), though we have not calculated them explicitly. 
This quantum conductance has to be compared to the classical conductance obtained from Eq. (4.4) by setting 
\r he \ 2 * T 2 /4: 

For both the Orthogonal and Symplcctic Ensembles the conductance is enhanced above the classical value, with a 
correction O(N). Note in particular that in the absence of a field the conductance increases as disorder is added to 
the junction. When the length of the junction is such that Ts ~ 1, the junction conductance reaches its maximum 
value of roughly g^s — NT/2. The classical value at this length is g^ s ~ NT 2 /2. This very large enhancement of 
the conductance is referred to as "reflectionless tunnelling." 

In practice it is the sensitivity of the reflectionless tunnelling effect to a magnetic field which is observed in ex- 
periments [4]. We have seen explicitly in the analysis of the quasi-ballistic junction that the phase conjugation of 
electrons and holes traversing time reversed paths is destroyed in a magnetic field, and only the O(N ) correction, 
due to interference involving identical paths, should remain. Thus, by applying a magnetic field, the conductance of 
the junction should be dramatically reduced to its classical value, within a correction of O(N) . This argument also 
indicates that if we break spin rotation symmetry, the qualitative behaviour in a magnetic field should be unchanged, 
and there should only be quantitative corrections at O(N ) due to the spin orbit interaction. To confirm this we have 
carried out a numerical simulation of a junction under the relevant conditions. 

The model used has been described in Section VII. The parameters used in the simulation are L x = L y = 6, L = 48 
and Ep = for which TV = 24 at 9 = and TV = 22 at 9 = tt/8. We have set W — 2 and estimate that the mean 
free path I ~ 12. The junction is in the diffusive (L z » I) quasi-ld (L x , L y < I, L x ,L y <C L z ) regime. In Figure 6 we 
present the data for the conductance versus magnetic field obtained in the simulations for the two values of spin orbit 
interaction parameter 9 indicated above. The magnetic field is given in terms of the flux <f> penetrating the normal 
part of the junction. At z=L there is a potential barrier of height Uq — 3. The barrier transparency is a function of 
Uq and also has a weak dependence on the the spin orbit parameter 9 and the magnetic flux <f>. The average value of 
T = 0.23 with a variation of about 10% over the parameter range used. At zero field and zero bias the conductance 
should be approximately: 

2e 2 NT e 2 

° NS ~ ~h T ** T <"» 

which is in reasonable agreement with the numerical data. On applying a magnetic field we expect to approach the 
classical value: 

G d 2e 2 N ^ e 2 

This is again in reasonable agreement with the numerical data. A strong reduction of the conductance is observed 
both when spin rotation symmetry is unbroken 9 = and when it is broken 9 = tt/8 by the spin orbit interaction. The 
critical field corresponds to a flux <f) — 20o penetrating the sample, in agreement with critical field expected for the 
suppression of the weak localisation effect in normal metals. The data confirm that, apart from possible corrections 
of O(N ), the reflectionless tunnelling effect is not sensitive to the spin orbit interaction. 



IX. CONCLUSIONS 



We have studied the effect of the breaking of spin rotation symmetry, as a result of an appreciable spin orbit 
interaction, in determining the quantum correction to the classical conductance of a disordered mesoscopic normal 
superconducting junction. 

The most striking result we have obtained concerns the NS junction in an applied magnetic field. We have found 
that even if time reversal symmetry is broken, the quantum correction is sensitive to the spin orbit interaction. A 
semiclassical analysis of a quasi-ballistic junction permitted the scattering processes, which interfere to produce this 
effect, to be identified. They involve scattering paths in which an electron and a hole traverse a loop in the normal 
part of the junction in the same sense, that which we have called here "identical paths" . 
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APPENDIX A: 



In this appendix we group together certain technical details concerning the scattering theory for model (3.1) which 
appears in the main text. We adopt the Dirac notation so that, for example, < xyza\ip >— ip(xyza). The particle 
current J p (z) through a cross section at z can be shown to be: 

Mz) = ±r<1>\D[z]\il>> (Al) 

where D is the operator: 

D[z] = \xyza X xyza\H e \xyz + lfi X xyz + l/i| 

— \xyz + lcr >< xyz + la\H e \xyzfi >< xyz[i\ (A2) 



By rewriting the Schrocdingcr equation in the contact in terms of a real space transfer matrix, and considering the 
orthogonality relations between left and right eigenvectors of this matrix, it is possible to show that the Bloch states 
satisfy the orthogonality relation: 



< m\D[z]\n >= if k m ^ k* 

In addition for real k the transverse wavefunctions may be normalised so that: 

< ±m|D[z]| ± n >= ±iS ntm 

Each channel then carrys a charge current of e/K. Here we use the notation that: 

< xyza\ ± m >= ip± m (xya) exp(±i£: m ) 



(A3) 



(A4) 



(A5) 



In the absence of a magnetic field the Hamiltonian H e will commute with the time reversal operator T . In this 
case if ip e is a solution to the scattering problem we may construct another solution Tift by operating on i[> e with T. 
Taking note of this we can derive the following condition on S e if the Hamiltonian is time reversal invariant: 



S 
-(ft 



-d 
d 



where d is a 2N x 2N symmetric unitary matrix with elements: 

d m ,n = -» ^ i)\ m (L x - x + 1, y, a) [[v zt ] (T , Ai exp(ik m ) - [v%^ exp(-ik n )] ip+ n (x, y, n) 



(A6) 



(A7) 



The condition on the scattering matrix when H e is time reversal invariant, can be simplified by an appropriate choice 
of basis. Since d is symmetric unitary it may be decomposed (not uniquely) as d = e T e with e unitary. If we do this 
and make a transformation to the new basis: 



1p' +n = ^ e n,m1p+r 



(A8) 



the matrix d becomes unity d = \2N- This transformation is legitimate since it leaves the current normalisation 
unchanged. The general relation between the electron and hole scattering matrices arrived at is, as indicated in the 
main text: 
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S h (+e,+B) = 



d T 
-d T 



S* e (-e,+B) 



-d* 
d* 



(A9) 



which again may be simplified by the transformation (A8). 

We now turn to the relation between the electron scattering matrix S e and the electron Green function G+. Let 
us suppose that H = H e + U with U a random potential. To each incoming state at the left | + n > we associate a 
scattering state \4>+n >'■ 



\<t>+ n >= I + n > +G+U\ + n> 



(A10) 



Supposing that the random potential U is nonzero only in the volume V between zo < z < z\ this can be transformed 
to: 



\<P+ n > =GtD[H e , Zl }\+n> -GtD[H e ,z -l]\+n> 
+ \xyza >< xyza\ + n > 



In z < zq + 1 the scattering state has the form: 

\4>+n >= | + n > 



m > 



(All) 



(A12) 



and in the region z > z\ — 1: 



\<t>+n >= y~][te]m,n\ + m > 



With the aid of (A3) and (A4) we now find 



[r, 
[t, 
[r 

K 



mn 
mn 
eimn 

< III I! 



< -m\D[z }GtD[z -l}\+n> 
+i < +m\D[ Zl - l]G+D[z - 1] | + n > 
<+m\D[z 1 -l]G+D[z 1 ]\-n> 
= +i < -m\D[z + l]G+D[ Zl }\ -n> 



(A13) 



(A14) 



The last two relations are obtained by considering an incoming wave | — n > from the right. In a similar way a relation 
between Sh and G^ can be derived. We note only the expression for the reflection matrix at the right: 



[r'hlmn = -i< -m\D[ Zl - l]G+D[ Zl }\ +n> 
and we now assume that the magnetic field is zero except in V. 



(A15) 



APPENDIX B: 



As mentioned in the text, the simplicity of (3.18) is related to the presence of p in Eq.(3.5). In most references dealing 
with random matrix approaches to conduction in disordered solids definitions, such as those of ref. [16], p does not 
appear in the equivalent of these equations. To avoid confusion we describe explicitly the transformation between the 
sets of definitions. 

In the absence of a magnetic field, the scattering matrix of Ref. [16] which we shall call S^ p (containing reflection 
and transmission matrices r^f p , t^ p , etc.) has the form: 

(Bl) 

and the property that: 

(B2) 

where the 2N x 2N antisymmetric unitary matrix Z satisfies (C3) and (C4). The coefficients a+, a' + introduced in 
Section IV are identical to those of Ref. [16], while the coefficients a_, a'_ are related by: 



r MP .IMP 

e e 
+MP „IMP 



(Sf p f = 



Z T 





S MP 


z 








Z T 





z 
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a MP = Za_ 
> MP = Za'_ 



a 



The relation between S matrices is then: 



S e = 



Z T 
1 



1 
z 



(B3) 



(B4) 



so that (B2) and (B4) give the relation (3.9). 

In the absence of spin-orbit coupling the two spin components would be independent and the reflection and trans- 
mission matrices would be block diagonal so that, in a suitable basis, S MP would have the structure: 



-.MP 



f e f e 

f e t' e 

t e r' e 

L f' 



(B5) 



Here, f e ,t e , ... are N x N matrices that correspond to the reflection and transmission matrices for either spin direction 
and have the usual properties, so that S^f p is symmetric. The S of this paper, obtained from the transformation 
(B4), would have the structure: 



-f e V 

f e 

t e 

o L -f' 



o 

t e ' 

K 





(B6) 



APPENDIX C: 



The parameterisations of S e appropriate for the various ensembles are discussed in Refs. [3,16,24,33,34]. For the 
Unitary II Ensemble S e may be parameterised: 

U2VT113 U2VR114 

where the u's are 2N x 2N unitary matrices, T = diag{T„; n = 1, . . . , 2iV} and R = I2N — T. (The parameters T n 
appearing here are the same as the eigenvalues of t e tl appearing in the conductance formulae of Section VI.) The 
unitary matrices are independently and uniformly distributed with respect to the invariant measure on the group 
U(2N) of 2N x 2N unitary matrices. 

In the Symplectic Ensemble, time reversal symmetry imposes further restrictions on the form of S e and the appro- 
priate parameterisation is: 



ui^/RZuJ —ui^/TZuT 
VTZuf u 2 \ r RZu T 



with u\ and u 2 uniformly distributed on U(2N) and Z is a fixed antisymmetric 2N x 2N matrix satisfying: 

[Z,T]=0 



(C2) 



(C3) 



-1 



2N 



(C4) 



If the spin orbit interaction is negligible the spin is conserved during the electron motion and the Hamiltonian and 
scattering matrices can be block diagonalised. Thus in the Unitary I Ensemble a basis can be found such that: 



-ui\/Ru 3 U1VTM4 

u\\/~Ruz u\\Ttu^ 

u 2 Vtu 3 -U2V&U4 

u 2 \/Tu 3 u 2 \/ r Ru4, 



(C5) 



18 



where u\, . . . , U4 are N x N unitary matrices uniformly distributed on the unitary group U(N) and f = diag{T„; n — 
1,...,N} and R= l N -f. 

Finally in the Orthogonal ensemble S e is parameterised: 



S K = 



-ui 
ui\/Rw[ 
u 2 VtuJ 



Rui u\ 







Ul 

— U 2 \fRuK 



u 2 VtuJ u 2 ^TRu% 



(C6) 



with ui and u 2 uniformly distributed on U(N). 

The distribution of the T„'s is given by the solution of a Fokker Planck equation, the precise form of which depends 
on the ensemble under consideration. The joint distribution, and in particular the correlations, of the T„'s depend 
on the symmetry of the ensemble under consideration. However, all we need to know for the present purpose is that: 
[35,36] 



27V 



]TT„) =2N(l- S ) + 0( S 2 ) 



(C7) 



\n=l 



independent of the ensemble. The only further information needed to carry out the required averages is that if u is 
an N x N unitary matrix uniformly distributed on U(N) then 



1 



Si.i'Si 



(C8) 



An important property of the local maximum entropy model is that the distribution of S e is isotropic, that is to say 
the matrices Mi, . . . , u 4 are distributed according to the invariant measure on the unitary group and independently of 
the parameters T\, . . . , T 2 n- 
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Figure captions: 

Figure 1: A schematic of the NS junction for which the scattering theory is developed in Section III. 

Figure 2: An example of a path which contributes to process p\ corresponding to an electron (solid line) traversing 
the normal part of the junction whose path is then retraced by the Andreev reflected hole (dashed line). 

Figure 3: An example of a "time reversed path" which contributes to process p$. The path of the electron moving 
from x a to Xb is retraced by the Andreev reflected hole as it moves from Xb to x a . Interference between this path and 
that illustrated in Figure 2 is insensitive to the spin orbit interaction and is suppressed in a magnetic field. 

Figure 4: An example of an "identical path" which contributes to the process p§. The electron moves around a loop 
and returns to x a . The Andreev reflected hole traverses the loop in the same direction. Interference between this 
path and that illustrated in Figure 2 is sensitive to the spin orbit interaction but not to a magnetic field. 

Figure 5: The conductance of a quasi-ld diffusive NS junction as a function of the spin orbit interaction parameter 9 
in the presence of an applied magnetic field ( flux of l/25^o per lattice cell). The barrier transparency T = 1. The 
NS conductance is sensitive to the breaking of spin rotation symmetry even though time reversal symmetry is broken, 
while the normal conductance is approximately constant as expected. 

Figure 6: The conductance of a quasi-l<i NS junction as a function of magnetic flux penetrating the normal part of 
the junction for zero 9 = and strong 9 = n/8 spin orbit interaction. The barrier transparency is low (r ~ 0.23) and 
the length of the junction is such that Ts ~ 1. We see that the qualitative features of the reflectionless tunnelling 
effect are insensitive to the spin orbit interaction. 



TABLE I. The various ensembles for which Gns is calculated. The abbreviation TRS means time reversal symmetry and 
SRS spin rotation symmetry. 



Ensemble 


TRS 


SRS 


Orthogonal 


yes 


yes 


Unitary I 


broken 


yes 


Unitary II 


broken 


broken 


Symplectic 


yes 


broken 



TABLE II. The quantum correction SGns = (2e 2 /h)SgNS for the quasi-ballistic NS junction for the ensembles listed in 
Table 1. 



Ensemble 


SgNS 


r = l 


r < l 


Orthogonal/Symplectic 


-2Nsf(T) 


-4Ns 


+iVsr 2 


Unitary I 


-2a/(r) 


-4s 


+ s r 2 


Unitary II 


+sf(T) 


+2s 


- s r 2 /2 



TABLE III. The quantum correction SGns = (2e 2 /h)SgNS for a diffusive (1 <C s <C N) NS junction without a Schottky 
barrier (r = 1). 



Ensemble 


SgNS 


Orthogonal 
Symplectic 
Unitary I 
Unitary II 


4 - 1 (-0.593) 
-1(4,- 1) (+0.297) 
-i (-0.333) 
+ | (+0.166) 
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